Transport Properties of One-Dimensional Hubbard Models 



O 

o 
o 

> 

O 



•4-* 



c 

O 

o 



> 
oo 

^- 

o 

oo 

On 
-i— > 



T3 

c 

o 
o 



X 



S. Kirchner 1 * 2 , H.G. Evertz 1 ' 3 , and W. Hanke 1 
1 Inst. f. Theoretische Physik, Univ. Wiirzburg, 97074 Wurzburg, Germany 
2 Inst. f. Theorie d. Kondensierten Materie, Univ. Karlsruhe, 76128 Karlsruhe, Germany, 
3 ISSP, Univ. of Tokyo, Roppongi 7-22-1, Minato-ku, Tokyo 106, Japan 

(April 15, 1998) 

We present results for the zero and finite temperature Drude weight D(T) and for the Meissner 
fraction of the attractive and the repulsive Hubbard model, as well as for the model with next nearest 
neighbor repulsion. They are based on Quantum Monte Carlo studies and on the Bethe ansatz. We 
show that the Drude weight is well defined as an extrapolation on the imaginary frequency axis, 
even for finite temperature. The temperature, filling, and system size dependence of D is obtained. 
We find counterexamples to a conjectured connection of dissipationless transport and integrability 
of lattice models. 



I. INTRODUCTION 

An ideal conductor is characterized at zero tempera- 
ture by a nonvanishing Drude weight D in the real part 
of the conductivity, 



Re{cr(w)} = D8{u) + a reg (iv) 



(1) 



as first introduced by Kohn in the context of the Mott 
transition p|. For finite temperature, the Drude weight 
D(T) can be introduced by a formal extension of eq. (|lj): 



Re{a(u;,T)} = D{T)S(u) + a reg {u,T) 



(2) 



A superconductor is characterized by an additional quan- 
tity probing the Meissner effect, the superfluid density 
||. There are similar transport quantities for spin de- 
grees of freedom . 

Recently there has been a lot of interest in transport 
properties at finite temperature, and especially in the 
Drude weight, since the High-T c materials exhibit an un- 
usual low-frequency behavior of the optical conductivity 
in the underdoped regime m. Based on analytical and 
numerical results, a very interesting possible connection 
between the integrability of a lattice-model and its finite- 
temperature Drude weight has been proposed, conjectur- 
ing that an integrable system is characterized by a 
finite Drude weight D(T > 0) ^ when D(T = 0) ^ 0, 
and remains an ideal insulator D(T > 0) = when 
D(T = Q) = 0, whereas a nonintegrable system should 
exhibit a vanishing Drude weight at T > 0. 

The results of our finite temperature Quantum Monte 
Carlo simulations do not confirm such a connection in the 
repulsive Hubbard model with and without next near- 
est neighbor interaction. We also show that the Drude 
weight can be extracted directly by an extrapolation of 
the current-current correlations on Matsubara frequen- 
cies, even at finite temperature, thus avoiding an analytic 
continuation to real frequencies. 

We investigate the one-dimensional Hubbard model on 
a ring of L sites threaded by a flux L$. The Hamiltonian 
is 



-*£< 



C i,a C i+l,a ( 



-ie <J>(a?i ,t) 



(3) 



where c, c* are annihilation and creation operators and 
the Peier Is phase 3>(i?j, u) = J s * +1 A(z, lu) dz in general is 
a function of position and frequency. We set h, c and the 
lattice spacing to unity, and we specify energies in units 
of t. We use periodic boundary conditions Q + i, ff = Ci >a . 

The Drude weight of the one-dimensional Hubbard 
model at zero temperature has been investigated in sev- 
eral papers, including studies of the scaling behavior of 
D at half-filling by Stafford, Millis, and Shastry ||, close 
to half-filling by Stafford and Millis (!(]] and by Fye et 
al. pT[ ] . For arbitrary filling, results were given by Schulz 
p"2[ , and by Fye et al. with the Lanczos method for 
small systems JTT| ]. Romer and Punnoose 113] computed 
both Drude weight and spin stiffness. Some related prop- 
erties of charge and spin currents at finite temperature 
were recently computed by Peres et al. [{lif in a per- 
turbation theory based on the Bethe ansatz. For the 
limit L — ► oo, expressions for the Drude weight based on 
the Bethe ansatz at finite temperature were very recently 
given by Fujimoto and Kawakami fjq| . 

In section [n] we discuss representations of the finite 
temperature Drude weight, and show that it can be ob- 
tained by an extrapolation purely in imaginary frequen- 
cies. In section III we briefly review the conjectured con- 
nection to integrability. In section IV we compute the 
filling dependence of D(T = 0) at zero temperature via 
the Bethe ansatz equations, and use these equations to 
give an approximation to D(T) at half-filling for low tem- 
perature and small system sizes. In section ^ we present 
a systematic Quantum Monte Carlo study of the finite 
temperature Drude weight in the Hubbard model, both 
repulsive and attractive. Results for a closely related 
property, the Meissner fraction, are presented in section 
N% and we study D(T) in the extended Hubbard model 
in section VII. 
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II. THE DRUDE WEIGHT 



Applying a Kramers-Kronig relation to eq. (0) results 



D(T) 



7r lim (u Im{er(aj, T)}). 

uj— >o 



Using linear response, the Drude weight is given by 



D(T) 



(k x ) - Re{A(0,w -> 0)} 



(4) 



(5) 



at $ = 0, where (k x ) denotes the average kinetic en- 
ergy per site and A(q,u>) is the current-current correla- 
tion function in frequency space ra] (see appendix A for 
details). 

There are two different ways to obtain the Drude weight 
from the imaginary frequency correlation function. One 
can perform an analytic continuation of the data to real 
frequency with the Maximum Entropy method to obtain 
ImA(w)/w, and then use the f-sum rule: 



ire 1 * 7r 



duj 



ImA(w) 



o+ 



(6) 



Alternatively, one can work entirely on the imagi- 
nary frequency axis: The analytical continuation of the 



current-current correlation function A, given in Eq. (A4), 



is valid in the continuous upper plane, including the 
imaginary axis at frequencies different from the Matsub- 
ara frequencies. One can therefore take the limit lj — > 
for A either along the real axis, or purely on the imagi- 
nary axis, even at finite temperature. The latter version 
eliminates the need for an analytic continuation (e.g. via 
Maximum Entropy) from data on the imaginary Mat- 
subara frequencies onto real frequencies. In the present 
paper we employed this procedure. We have verified that 
it produces results compatible with those using the f-sum 
rule, but with smaller errors. 

The generalization of D to finite temperatures can also 
be achieved by defining H 



e -f3E n d 2 E n (<S>) 



d $ 2 



$=0,ij=0 



(7) 



where Z = tr e @ H is the partition function. In the limit 
T — > 0, this immediately reduces to Kohn's Drude weight 



. ird 2 E (<S>) 



(8) 



The equivalence of eq. (|^) and eq. (^) within perturbation 
theory is shown in appendix A. 



III. CONNECTION TO INTEGRABILITY 

A model is usually called integrable if the energy eigen- 
values are distributed according to a Poisson distribution. 



For a non-integrable model, the eigenvalues follow a GOE 
distribution (Gaussian Orthogonal Ensemble) JFtJ . It can 
be shown for lattice models which are solvable by the 
Bethe ansatz that the eigenvalues are Poisson distributed 
]l8[ . Hence the Hubbard model eq. (||) is integrable. In 



section VII we will study the extended Hubbard-modcl 
with nearest neighbor repulsion, which is non-integrable 

& 

Based on analytical and numerical results, Zotos et 
al. 0-^1 have conjectured a very interesting connec- 
tion between the integrability of a lattice model and its 
finite-temperature Drude weight, stating that for a one- 
dimensional model the finite temperature Drude weight 
in the thermodynamic limit is 

(1) nonzero for an integrable system when it is nonzero 
at T = 0, 

(2) zero for an integrable system when it is zero at 
T = 0, and 

(3) zero for a non-integrable system. 

Conjecture (2) is different from the original suggestion 
jl| of a direct equivalence between integrability and fi- 
nite Drude weight at T > in one-dimensional systems. 
It was motivated by an explicit example in which the au- 
thors computed the Drude weight for a one-dimensional 
model of interacting fermions in the insulating regime 
with a "Mott-Hubbard" gap, and found a vanishing 
Drude value in the insulating regime even at high tem- 
peratures. 

Note that there is no rigorous proof of a connection be- 
tween integrability and a finite Drude weight. However, 
subject to some assumptions, Zotos et al. have shown in 
their most recent paper || that the Drude weight is finite 
whenever there exists an operator A such that [A, H] = 
and 



lim 



L (A 2 ) 



?0. 



(9) 



Thus, a weaker condition than integrability might suffice 
to make the Drude weight finite. 



IV. ZERO AND LOW TEMPERATURE: BETHE 
ANSATZ 

We employ the Bethe ansatz equations for the Hub- 
bard model, first obtained by Lieb und Wu in 1968 [pof , 
to calculate D(T) and its filling dependence exactly at 
zero temperature, and to provide an approximation for 
D(T) at low finite temperature. As was shown by Shas- 
try and Sutherland 0] , the Hubbard model with twisted 
periodic boundary conditions 



*(..., r + Lx, ...) = exp[i0] *(..., r, . 



(10) 



2 



can be solved with an appropriate ansatz for the wave 
function. Here <f> = eL $ is the overall phase aquired 
along the chain of length L |l6| . The Bethe ansatz equa- 
tions derived by Shastry and Sutherland are: 



M 



Lk n = 27r/ n + + 2 arctan[4(Aj- - smk n )/U] 
3=1 

(11) 



N 



2 arctan 

71=1 



[4(Aj- sin £;„)/[/] 



(12) 



A I 



= 2wJj + 2 ^ arctan[2(Aj - K)/U], 

where N is the total number of electrons, M is the number 
of electrons with spin down, the A,,- are parameters asso- 
ciated with the spin dynamics, and the quantum numbers 
I n and J m characterize charge and spin excitations. The 
quantum numbers have to be chosen such that [10 f 



/„ is integer/half-integer if M is even/odd, 

J m is integer/half-integer if N — M is odd/even, 



all I n (n — 1, .., N) (all J m (m = 1, 
from each other, and 

\In\ < (L-l)/2 

I J™ I < (N-M)/2. 



, M)) are different 
(13) 



The energy of a state with given quantum numbers is 
given by pi 



Following ref. 
from eq. (]?]) as 



D = 27re 2 L^ < cosfcn 



E = -2 ^cos(fc„) . (14) 

n 

we then calculate the Drude weight 



sin kr, 



= 

(15) 



A. Zero Temperature 

At half-filling and zero temperature the asymptotic In- 
dependence of the Drudc-weight is 

D(L;T = Q) = -(-l) L / 2 C(U)L l / 2 exp[-L/((U)}, (16) 

(L even) as was shown by Stafford and Millis |u|. The 
correlation length C is defined by eq. Jl6|). For L — > oo the 
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FIG. 1. Scaling dependence of D(T = 0) at U = 4 and 
half-filling. 
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FIG. 2. The Drude weight at U = 4 and 

quarter filling i/ = | scales towards a finite value 
D(L = oo)/(2e 2 ) = 1.29. The inset shows the scaling depen- 
dence (D(L = oo) - D(L))/(2e 2 ) = L _1 ' 12(1) exp[-L ■ 0.04(1)] 
on a logarithmic scale (at L mod 4 = 0). 



Drude weight in eq. ([L6|) vanishes, indicating the finite 
charge gap at half-filling. We show results for half-filling 
in fig. For L mod 4 = the ring becomes paramag- 
netic. This has already been discussed in refs. || and [[ll] . 
We find that the effective correlation length £ ss 4.25 at 
the system sizes of fig. [j] is slightly larger than the value 
C = 4.06 at U = 4 found numerically in |J. 

Away from half-filling the Hubbard model describes an 
ideal metal, with a finite Drude weight at L — > oo. We 
show the system size dependence at quarter filling in fig. 
||. We find that the L-dependence is very similar to the 
half- filled case, shifted by D(L = oo). Even though eq. 
( |l6| ) was only derived for the half-filled case, it fits the 
approach to D(L — oo) very well, albeit with a different 
exponent for L and a very large value of £, as shown in 
the inset of fig. ^. 
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1.5 



1.0 



Inverse Temperature 


a 


c 


P = 60 


16.65 ±0.16 


0.2400 ± 0.0003 


P = 80 


11.90 ± 0.11 


0.2400 ± 0.0003 


p = 100 


7.140 ±0.068 


0.2400 ± 0.0003 


P = 120 


2.379 ± 0.023 


0.2400 ± 0.0003 



TABLE I. Parameters a and c obtained by fitting the 
curves in fig. H to eq. (|l8h. 
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FIG. 3. The filling dependence of D at U = 4 and T = 0, 
extrapolated to infinite system size. 



We determined the complete filling dependence of the 
Drude weight by solving the Bethe ansatz equations nu- 
merically for system sizes up to L = 80, and extrapolat- 
ing to infinite system size using eq. ( |l6| ) . The results are 
shown in fig. |[ As expected, the Drude weight vanishes 
for v = (no electrons) and v = 1 (L electrons). For 
small system sizes up to L — 12 the filling dependence 
has already been reported by Fye et al. using Lanczos 
techniques jfl||. Solving the Bethe ansatz equations en- 
ables us to access much larger systems (up to L = 80) 
and thus to obtain a better extrapolation to the thermo- 
dynamic limit. The Drude weight for fixed system size 
L = 100 has also been obtained by Romer and Punnoose 
[Of, and our results agree. 



B. Partial Drude weight at finite temperature 

We now turn to finite temperatures and half-filling. 
Whereas ground state properties depending only on the 
energy can be calculated easily using the Bethe ansatz 
equation eq. (|i"T|), dynamical properties are much harder 
to determine via the usual Bethe ansatz. One recent 
method by Klumper et al. reduces the problem to finding 
the solution of two coupled nonlinear integral equations 
p3| . Here we will be content with a much simpler pro- 
cedure to calculate the spin triplet contribution to the 
finite temperature Drude weight at half-filling. We ne- 
glect excitations across the charge gap. As we will see, 
this approximation is valid for small values of LT/U, i.e. 
for low temperatures and not too large systems. 

Low energy excitations are described by variations of 
the quantum numbers around the ground state configu- 
ration. We choose N mod 4 = 2 whence the ground state 
is a singlet [[l0| . As excitations around this state we con- 
sider spin triplet states, with M — N/2 — 1, that is, we 



consider the 5=1 and S z = — 1 excitations. To include 
singlet excitations one has to refer to the general Bethe 
ansatz equations introduced by Takahashi p4[ , with a 
phase included as in ref. jl4| . The generalized equations 
describing all four branches of singlet and triplet exci- 
tations (the afore mentioned ones and the S — 1 and 
S z = excitations) are given in J2f]]. 
Due to the charge gap at half-filling, there are no low 
lying charge exitations. The charge quantum numbers 
I n therefore remain unchanged from the ground state. 
For the spin quantum numbers we get from eq. ( |i~3| ) with 
M = N/2-1 that 



\Jj\ < 



N 

T 



(17) 



and the Jj are now half- integer. Thus there are N/2 ± 1 
possible values for the M = N/2 — 1 spin quantum num- 
bers. All possibilities to distribute the J,'s according to 
eq. ([l7]) that differ from the ground state constellation 
describe excited states. These are the lowest possible 
excitations. In addition, there is a singlet excitation, de- 
generate in energy with the three triplet excitations [^2) . 
This singlet excitation also contains two holes in the spin 
quantum numbers compared with the ground state con- 
stellation. It is described by complex quantum numbers. 
Due to the degeneracy, its contribution just multiplies 
the partition function by a constant factor. 

The quantity D tri P let {/3, L) calculated via eq. (|l|) us- 
ing only the triplet spin excitations is of course only part 
of the Drude weight, but as we will see it should be a good 
approximation for small temperatures and small system 
sizes. 

For the temperature and system size dependence of 
■Qtnpiet^p^ jjj ^ ^ moc } 4 — 2) we find again a relation 
similar to the T = result of Stafford and Millis |l(J , eq. 



D triplet (/3,L) 



2e 2 



= a(P)L b exp[— c • L], 



(18) 



where only a((3) is a function of the temperature. The 
results for ]j tn P let are shown in fig. ^, and the fitted 
parameters a and c in table |. 

We conclude that in the thermodynamic limit, L — » 00, 
triplet excitations do not contribute to the Drude weight. 
This also remains true for higher spin excitations [p6| . 
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V. FINITE TEMPERATURE: QMC RESULTS 
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FIG. 4. Scaling behavior of D triplet (T) as a function of L 
for finite temperatures and U = 4, v = 1. 



For the /3-range studied, the exponent b = 0.336±0.005 
and the length scale c _1 = 4.17 ± 0.01 are independent 
of temperature. They differ somewhat in value from the 
T = result for the complete Drude weight D (T = 0) 
(see eq. © and fig. Q). The amplitude a(/3) of D tri P let 
is precisely linear in (3 for the /3-range studied. 

a(f3) = (30.911 ± 0.007) - (0.23771 ± 0.00009 ) • (3. (19) 

This is surprising, since it resembles the high tempera- 
ture behavior generally expected of the complete Drude 
weight §. 

We found that the interaction dependence of D triplet 
for a given temperature and system size can be approxi- 
mated by 



D tnplet (U) ~ U- d • exp[-/ • U] 



(20) 



where the parameters d and / depend on system size. 

Recently we received a paper by Fujimoto et al. |l5| ] 
with exact expressions for the Drude weight at finite 
temperature in the limit L — > oo, based on the Bethe 
ansatz. The Drude weight is finite at finite temperature, 
in agreement with our results. The authors compute ex- 
plicitely the leading contributions to ~D(L — oo,T) at 
small temperatures, which at half filling are, as expected, 
proportional to exp[— AmhP], where Amh is the Mott- 
Hubbard gap. 

Since the quantity D t7 " lplet ((3 , L) vanishes exponen- 
tially as I — * oo, the spin contributions to the Drude 
weight considered here should dominate over the charge 
excitations for mesoscopic systems with (3 ^ 1 (so that 
higher excitations are suppressed) and 

cL < PU (21) 

i.e. for low temperatures T and small values of LT/U. 



In this section we present our results for the Drude 
weight at large finite temperature and compare them 
with the prediction for the temperature and system size 
dependence obtained in the previous section. 

The simulations were carried out using the grand 
canonical Quantum Monte-Carlo method |27|] . The cal- 
culation yields the current-current correlation function at 
discrete imaginary times r^. Performing a fourier trans- 
formation results in the current-current correlation func- 
tion at the bosonic Matsubara frequencies iuj n . The an- 
alytic continuation of this function onto real frequencies 
at zero momentum gives A(0,o;) of eq. (||). However, 
the analytic continuation is a numerically ill-conditioned 
problem, which we can avoid. As pointed out in section 
P, one can take the limit ui — > directly along the imag- 
inary axis. We therefore calculate D by fitting A on the 
Matsubara frequencies. This is still an ill-conditioned 
problem when only information about a finite number of 
Matsubara frequencies is available. It is therefore impor- 
tant to use a fitting function with the proper analytical 
form. 



From the analytic continuation of eq. (A4) it can be 
seen that A(0, iuji) is a well-behaved function on the imag- 
inary axis, namely a sum of Lorentz curves: 



A(0,iwj) 



3 



3 , i 2 



We approximate it by a finite series 

n r 

Afit (0,iwj) = - 



h 2 



d 2 



+ P 



(22) 



(23) 



-fit. 



and determine the constants a,b,c,d,e,f by a \ 
This ansatz fits the data very well. Fig. |5| shows an ex- 
ample of our fit. We find that the third term in eq. (|23|), 
e/(ojf + f 2 ), and often also the second term, do not con- 
tribute significantly, indicating that eq. ( ^3| ) is a good 
approximation. However, it should be noted that the 
constants a, c and e are not identical to the differences of 
the energy eigenvalues of H appearing in A(0, iuji), since 
they result from a fit to a truncated series. We have ver- 
ified that our procedure to calculate the Drude weight 
produces results compatible with those using the f-sum 
rule eq. ([}]), but with smaller errors. 

We used this procedure to determine D(T) from finite 
temperature Quantum Monte Carlo runs. We typically 
collected about 180000 Monte Carlo sweeps for each data 
point, with a discrete time step of At = 1/32 for (3 — 3 
and At = 1/10 for (3 = 8. We obtained the following 
results. 

Fig. [6] shows the Drude weight for repulsive interac- 
tion U = 4 and inverse temperature (3 — 8. The system 
size dependence is compatible to that of D triplet (L) , eq. 
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FIG. 5. Fit of A(0,iu>„), at the nonzero Matsubara fre- 
quencies, for L = 32, U = 4, /i = 2.2 and /? = 10, plotted vs 
n. The cross marks the extrapolated value at u = 0. 
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FIG. 6. The Drude weight at (3 = 8, U = 4, v = 1 is 
nonzero in the thermodynamic limit L — > oo. The form of 
the interpolating function is motivated by the T = Bethe 
ansatz results. 

( |l8| ) , as expected for temperatures small compared to the 
charge gap. 

D(L)/(ire 2 ) = -A ■ L~ B exp[-C • L] + d. (24) 

We determined the parameters by a x 2_n t- For the case 
shown in fig. | the result is A = 10.12 ± 1.14, B = 
0.75 ± 0.88, C = 0.22 ± 0.09, and d = 0.020 ± 0.003. 

Thus, in the thermodynamic limit, we obtain a fi- 
nite Drude weight D{j3 = 8, v = 1, U = 4)/(7re 2 ) = 
0.02 ± 0.003. The system size behavior of the full Drude 
weight D(L, = 8) (eq. (§4|)) appears to be similar to 
that of the partial Drude weight D t " plet (L , f3) obtained 
by Bethe ansatz in eq. (|l^), within the large error of B, 
even though D(L,f3 — 8) contains charge exitations, as 
well as higher spin and charge-spin-excitations. 

Towards higher temperatures, the Drude weight in- 
creases rapidly. Results at /3 = 3 are shown in fig. |t]. 
The approach to the asymptotic value is very fast now, 
and there is no longer a clear exponential behavior. 

We checked that for lattice sizes L = An + 2 the finite 
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FIG. 7. At large temperature (/3 = 3, U = 4, v = 1) the 
Drude weight is large and converges rapidly with system size. 
The Meissner fraction k (see section WW) converges to zero 
quickly. 
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FIG. 8. The Drude weight in the repulsive model 
(U = 4,f — 1) increases rapidly with temperature. 

value of the Drude weight remains the same within errors. 
One might be concerned that at temperature (3 = 3 and 
lattice size L — 40 the effects of the finite size gap could 
still be important. We therefore verified that at very 
high temperature ((3 = 0.5, L = 20) the Drude weight 
remains finite. In fact, its value increases further, to 
about D/(ire 2 ) = 0.23(1). In fig. | our results for the 
Drude weight at various temperature in the half-filled 
case are shown. 

Away from half-filling, the finite temperature Drude 
weight remains finite. Results at quarter filling are shown 
in fig. ^. The asymptotic value is already reached for very 
small system sizes. Note that away from half-filling, there 
was some fluctuation of the actual filling for different L 
in the grand canonical QMC, which is not reflected in the 
error bars. This can cause the outlying point at L = 28 
in fig. |. 

For the attractive Hubbard model at half-filling the situ- 
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FIG. 9. At quarter filling v ^ 0.500 (/3 = 3, U = 4) 
the Drude weight is finite. Small variations in filling of the 
grand canonical simulations may be responsible for the large 
deviation of the data at L — 28. 
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FIG. 11. The Drude weight in the attractive model at 
lower temperature: /3 = 10, U = —4, v = 1 shows strong 
finite size dependence. 

VI. THE MEISSNER FRACTION 
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FIG. 10. Attractive Hubbard model at high temperature: 
(3 = 3, U = —4, v = 1. The Drude weight (upper points) 
quickly converges to a large finite value, whereas the Meiss- 
ner fraction k vanishes (lower points; see section |v| ). For 
comparison: (k x ) ~ —0.890, with little size dependence. 



ation is similar to the repulsive case. At low temperatures 
the system size dependence is slower, as shown in fig. 
and the Drude weight has not yet reached its asymptotic 
value in our calculations. At high temperature, shown 



in fig. 10, the Drude weight converges quickly to a large 
value. 

We see that in the half filled Hubbard model at finite 
temperature the Drude weight is non-zero in the thermo- 
dynamic limit L — > oo, both for the repulsive and for the 
attractive case. Since half-filling is the insulating case of 
the Hubbard-model at zero temperature, this result is in 



disagreement with conjecture (2) (section HI) 



A property closely related to the Drude weight is the 
Meissner fraction k lEHl defined as 



nd 2 F 



K — — 



Ld<P 2 



*=o,<f-^o 



(25) 



where F is the free energy. It reduces to Kohn's Drude 
weight for T = 0. 

Evaluation of eq. (J25|) results in 



(k x ) - A(<f-> 0,iw n = 0). 



(26) 



Note that A(q — > 0, iu> n = 0) at the zeroeth Matsubara 
frequency, given in eq. (A5), is in general not equal to 
the analytically continued correlation function at u> — > 0, 



eq. (A6), so that at finite temperature k differs from the 
Drude weight D by the contributions from degenerate 
states 



D{T) 



lim 

9->o L 



E Pn\(n\f(q)\m)\ 2 



(27) 



E m =E„ 



In the thermodynamic limit L — > oo and in a transverse 
vector field, k measures the superfluid density |0,lS,E8[. 
In our one-dimensional case it will thus become zero. For 
a finite system, the Meissner fraction can be finite. 

We show results for the Meissner fraction on finite sys- 
tems in figures [?],[l(],|l^ and [l3|. For L — > oo, k van- 
ishes. We find a finite size behavior similar to that of the 
Drude weight. As for the Drude weight at T — (see 
fig. ^), k is positive when L = An + 2 and negative when 
L = An. At high temperature = 3, k converges to zero 
very quickly both for the repulsive and for the attractive 
model (figures |^ and |l0|). At low temperature the ap- 
proach to zero is exponential in system size, but much 
slower, as shown in figures ^| and For very small 
systems, k is of similar magnitude as the Drude weight. 
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FIG. 12. Meissner fraction in the repulsive model at low 
temperature: f) — 8 and [7 = 4, (k x ) w -0.955, v = 1. The 
two values at 1/ = 10 and L = 26 differ in sign, since they 
belong to system sizes where L — 4 • n + 2. 




FIG. 13. Meissner fraction in the attractive model 
at low temperature, (3 = 10, U = —4, v = 1 
exp[-(0.1 ± 0.002) ■ L]. The kinetic energy is (fc^) ? 
for all system sizes. 



-0.96 



Even though the difference D — k, eq. Q27j), contains a fac- 
tor of /?, this difference turns out to be smaller at larger 
/? in the temperature range studied, both for small and 
for large systems, contrary to previous expectations p8[. 



VII. EXTENDED HUBBARD MODEL 

In section |V| we showed that for the integrable Hubbard 
model eq. (|3|) the finite temperature Drude weight both 
at half-filling and away from it is nonzero even in the 
thermodynamic limit. This is in contrast to the T = 
case, where the Drude weight vanishes at half-filling due 
to the charge gap. 

Let us now examine a nonintegrable case, provided by 
the extended Hubbard model with nearest neighbor in- 
teraction lOl 



0.6 
0.5 
0.4 



S °- 3 
£5 



0.2 
0.1 
0.0 



15 20 25 

system size L 



FIG. 14. The Drude weight in the extended Hubbard 
model at at U = 4, V = 1/2, = 3 and v = 0.8978 ±0.0002 is 
nearly constant. The error bars do not take into account that 
the filling has some variation with L in our grand canonical 
simulations. 



We have measured the Drude weight and the Meissner 
fraction at finite temperature in this model. Results for 
a weak extended interaction of V = 1/2 and away from 
half-filling are shown in fig. 
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H P = H 



H 



v z^( n n + n »i)( n »+it + m+u)- ( 28 ) 



The Meissner fraction 
drops to zero quickly, similar to the behavior in fig. [t]. 
The Drude weight has a large value already at small sys- 
tem sizes. It shows a very slow falloff, consistent with 
D ~ e -(o.ooo8i±o.oooo4)^ Thig faUoff may entirely be 

caused by a very small systematic increase in the fill- 
ing factor measured in the grand canonical simulations. 
If the exponential falloff is real, then the Drude weight 
would drop to small values only at extremely large system 
sizes, but would be zero in the thermodynamic limit, in 
agreement with the conjecture of Zotos et al.. However, 
from our data it appears more likely that the thermody- 
namic limit of D is finite, which would be in disagreement 
with conjecture (3). This would imply that, even though 
the model is not integrable, there is an operator satisfy- 
ing eq. (^) in the extended Hubbard model. A similar 
conclusion has been drawn for a related model recently 
©■ 

We have also measured the Drude weight at (3 = 5, U = 
4 and larger repulsion V = 1.6, on a large system L = 60. 
Both at quarter filling v = 0.5 and at v = 0.633 we find 
a finite value D/(ne 2 ) — 0.21(2), even though the model 
is not integrable. 



VIII. CONCLUSIONS 

We have presented results for the Drude weight D(T), 
on mesoscopic systems and extrapolated to the thermo- 
dynamic limit, for zero and finite temperature in several 
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versions of the ID Hubbard model, including the depen- 
dence on system size, filling, and temperature. 

At zero temperature we provided the full filling de- 
pendence of -D, by solving the Bethe ansatz equations 
on large systems. For small finite temperatures we com- 
puted the dominant contribution to D(T, L) on small sys- 
tems from the Bethe ansatz. At larger finite temperature 
we used Quantum Monte Carlo. We showed that the 
Drude weight can be obtained by an extrapolation of the 
current-current correlation function purely in imaginary 
frequencies. We found a non-vanishing Drude weight at 
finite temperature for all cases considered, the repulsive 
Hubbard model both at and away from half-filling, the 
attractive model at half filling, and the extended Hub- 
bard model away from half filling. The Drude weight 
quickly grows with temperature for the half filled Hub- 
bard model. 

Our results for the integrable half-filled Hubbard 
model do not confirm conjecture (2) (section III) on the 
connection between integrability and finite temperature 
Drude weight, and we find that conjecture (3) would dis- 
agree with the (likely) non-zero value of the Drude weight 
in the nonintegrable extended Hubbard-model. 
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APPENDIX A: SPECTRAL DECOMPOSITION 
OF THE DRUDE WEIGHT 



The Drude weight is given by 



D(T) 



ire 



= - (k x ) - Rc{A(0,cj^0)} 



(Al) 



at $ = 0. Here (k x ) is the average kinetic energy per site 
and A(g, u>) is the current-current correlation function in 
frequency space ||. It can be obtained from analytic 
continuation of 



A((f, iuj ri 







(f(q,r)j p (-qM dr (A2) 



where (3 — 1/T, iw n — 2irin/P are the Matsubara fre- 
quencies, q is the momentum of the applied external vec- 
tor potential, and ej p {q) is the Fourier transform of the 
paramagnetic current density, 



it 



E 

k = 1, ..,L 
<r=T4 



-lqx k 



t 

fe+1, 



• Cfe,, 



(A3) 

In an eigenbasis of the Hamiltonian eq. (|^), the current- 
current correlation function at nonzero frequency is given 
by 



A(q, iwi ^ 0) = 



E ^—v^ — ^\{n\f(q)\m)\\ 



L iui + E n — E>, 



ill 



(A4) 



and exactly at zero Matsubara frequency it is 



A(g,m> z = 0) 



E 



2_ 

L 

n, m 
E n j4 E„ 

! E 

n, m 
E n = E n 



E m -E„ 



\{n\f{q)\m)Y 



\{n\f{q)\m)f 



(A5) 



Here {n\j p \m) are matrix elements of the current opera- 
tor in the eigenbasis and P n = exp(—f3E n ) denotes the 
Boltzman factor for the nth eigenvalue of the Hamilto- 
nian. 

The current-current correlation function is well defined 
|]30f in the upper complex plane by specifying its val- 
ues on an infinite set of finite Matsubara frequencies, eq. 
( A4 ) . The analytic continuation can easily be performed 
by taking iu>i — > lu + id in eq. (A4). Performing the zero 
frequency limit u> — > then yields 



Re(A(g,o;^Q)) 



E 



Pn 



E m jtE n 



\{n\f{q)\m)\ 2 , 
(A6) 



which is identical to the first term in eq. (|A5|). Thus 
L 



E m 'jiE n 



E 



|(n|/(0)|m)| 2 



(A7) 



Using second order perturbation theory and assuming 
that $ removes all degeneracies, the expression takes the 
same form as eq. (H) , so that then the two definitions eq. 
(||) and eq. @) agree. 



The second term in eq. ( A5) is a nonanalytical part of the 
thermal Greens function which does not contribute 
to the Drude weight, but instead to the Meissner fraction, 
discussed in section VI. 
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